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An optical four-level atomic discrete system through optical induction is proposed. A theoretical 

scheme to produce nonclassical lattice solitons (NLS) in the system is presented with the use of 

the effects of enhanced self-phase modulation and the giant kerr effect in the electromagnetically 

C^ ' induced transparency. The power density and the photon flux can be tuned to a very low level by 

^—^ I the controlling field and the soliton can propagate with very slow group velocity. By changing the 

c7^ . sign of the detuning Ai , both in-phase and tt out-of-phase NLSs can be produced in this system. 

(N; 

^ \ Optical fields propagating in couple waveguide arrays exhibit novel phenomena [l|. The light field shows great 

^~~i I functionality in such an optical periodic system which are impossible appeared in the bulk. Recently, the behavior of 

f~^ nonclassical light field in couple waveguide arrays becomes an interesting issue [2, Q • Nonclassical light is always a 

light field with single or few photons. In quantum information and quantum communication, photons are widely used 

' jy. ' as the information carrier to transfer signals. Generally, the number of photons are very few |[^]. Therefore, such a 

Q I optical periodic system provides an avenue to inventing new optical devices for quantum information processing. To 

'^ ' avoid signals degradation during propagation, it is highly desirable to propagate the optical wave as a kind of soliton. 

Qh Discrete or lattice soliton (LS) [5, &\ is a soliton mode in such a periodic waveguide arrays. It shows great potential 

O in photonic network communications [7|. However, in the nonclassical light limit, the intensity of the soliton must 

c/3 be very low (i.e. at least few photons per soliton cross section per nanosecond). Traditional discrete systems such as 

kerr nonlinearity media [8|, liquid crystals [9] photorefractive crystals [lOJ are very difficult to produce such kind of 



O 



^ ■ soliton because a required large nonlinear refractive index is always associated with intense light field. Therefore, it 



Oh, 



is important to realize a new optical periodic system that can produce a good controllability NLS in guiding these 
signals. 

In this work, we will apply the technique of electromagnetically induce transparency (EIT) to overcome this difficulty 
in the atomic system. EIT media, which can be applied to inducing dissipation-free strong photon-photon interactions, 
can offer far larger third-order Kerr susceptibility and far smaller linear susceptibility than conventional nonlinear 
media [11]. Recently, ultraweak [12| and ultraslow [IS^ spatial and temporal solitons are realized by utilizing EIT in 
uniform atoms system, respectively. This opens the possibility of archiving NLS in atomic discrete systems. 

In EIT medium, a dissipation-free optical discrete system can be generated by using optical induction technique. 
The energy diagram of the scheme is shown in FIG. 1 (a). It is similar to the traditional N-type system: the parties 
of |1), |2) are same and different from |3), |4). A weak probe wave Ep with Rabi-frequency Vlp = p^iEp/H is acting 
i"-^ on a resonant transition |1) -^ |3) with a single photon detuning Ai. A running- wave field with a Rabi-frequency 
f~^ I flc is driving the atomic transition |2) -^ |3) with a detuning Ac = Ai, in other works, the two-photon detuning 
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FIG. 1: (a) Energy diagram of the system, (b) Geometric configuration of the lights. Lattice-forming beams fis and nonclassical 
light beam Ep are co-propagating. The coupling light Sic is a p-wave,which is propagating transversely along the atoms 
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(5 = Ai — Ac — 0. An optical induction field with a Rabi-frequency ils is inducing the transition |2) —^ |4) with a 
detuning A2. According to the previous studies [ll|, when the probe exactly at the two-photon resonance (<5 = 0) 
and the atoms are all in the ground state, the absorption to the probe is vanished and independent to the single- 
photon detuning. Meanwhile, the existence of the detuning Ai and A2 makes the probe experience an enhanced SPM 
produced by itself [l_4j and XPM by fls ^5|], respectively. 

We assume that rip^fls ^ ^c, which results in the action of flp and fls can be treat as a perturbation. The 
interaction Hamiltonian between the light and atoms can be write as: H = Hq -t- Hi, where 



i^o = E^^I*)(*|-9["ce-'"^*|3)(2|-fi/.c], 



(1) 



4=1 



is the Hamiltonian operator without perturbation and 



Hi 



1 



[npe-"^'''\3){l\ + ns\A){2\ + H.c]^ 
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is the perturbation's Hamiltonian. Therefore, the evolution equation of density matrix reads: 
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Where n designates the n-th step of the density matrix. The matrix elements (Jlr|j) = ^iSij, where 7^ is the decay 
rate designating the population damping from the energy level \i). Assuming that all the atoms are populating at 

the ground state, in other words, Pi-{ — 1 and p\, = (i or j 7^ 1), and then the steady state solution of any step 
of the density matrix can be deduced by Eq. (3). Here we only deduce to the 3-step, the higher steps of the density 
matrix are ignored for their influences are tiny to the system. In the deduction, we assume that 72 w 71 w and 
73 « 74 = 7. Furthermore, to avoid the absorption from the atoms, we let Ai, A2 S> 7, fie- From Eq. (l)-(3) and 
the initial condition of p^^\ one can obtain that: 
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where P21 



P12 ~ — i.Lp\lf-i i\\l(J 

And then obtain that: 
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where 0)42 — A2 — 1742 and 0)43 — A2 — Ai — ^743. The paraxial steady-state propagation equations of Ep and Es in 
the slowly varying envelopes read as: 
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FIG. 2: (a)In-phase LS with k = -2.5 (Ai = -IOO7) and Vb 
K = 2.5 (Ai = IOO7), Va = 0.75 and P = 0.0106. 



0.5 and P = / \Vfdi = 0.023. (b) tt out-of-phase LS with 



where iV is the density of the atoms. According to Eq. 
be treat as hnear by the condition of Ai, A2 ^ 7j, ^c- 
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(7), if we let |Ai| = IA2I, the propagation dynamics of Es can 
The reason is: (l)when Ai = A2, it is obvious obtained that 

pfj = 0; (2) when Ai = — A2, P^^-j I P^i '^ l^cP/^i ■^ 1; which resuhs in that the refractive index modulation from 
Ep can be neglected. Under this circumstances, Ep propagates nonlinearly with the modulation from Es, while Eq 
propagates nearly linearly. This relationship is the key element to form a discrete system via optical induction [5|. 
Therefore, we can apply \Es\'^ to form the pattern of the lattice waveguide arrays (See in FIG. 1 (b)). Substituting 

(3) 
P3I 



in Eq. (8), we have: 



where V^ — ( 
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i-^U = -^ViC/ + Vi^, ri)U + k\U\^U, 
with C, = zkci, ^ = xkci, rj = j/fco! U = Q,p/Q,c is the dimensionless field, and. 
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And then, Ep experiences two kinds of refractive index modifications: (1) a periodical index changes induced by f2s 
(via XPM); (2) an index changes induced by Ep itself (via SPM). Here, fts could be viewed as a linear potential 
added on Ep. The sign of detuning Ai defines that Ep experiencing self-focusing (Ai < 0) or self-defocusing Ai > 0; 
while the sign of detuning A2 defines that the energy of the soliton field are located on a lattice waveguide site (when 
A2 > 0) or between lattice sites (when A2 < 0). 

In the ID case, choose the lattices- forming wave to be \ils\^ — l^^sP sin^ (Tr^/d) , and the periodical potential in Eq. 
(11) can be written as V^ = V^ sin^(7r^/d), where Vq — N\p3i\'^\V,^\'^ /2eohA2\flc\^ ■ Assuming the probe experiences 
self-focusing (Ai < 0), the bound state solution of Eq. (10) can be obtained by the imaginary-time propagation (ITP) 
methods. The parameters are chosen referring to the Rb atoms as: p3i — 3.0 x 10"^^ D, 7 w 6 MHz and Xp ~ 795 nm. 
For convenience, other parameters are chosen as: d = 10, |i7g|^/|J7(7p = 0.1 with ^c — 27, Ai = A2 = — IOO7 and 
iV = 3 X 10^^ cm~3^ which are tunable in practice. Then the normalized parameters a and Vq are obtained: 
and Vo = 0.5. In the simulations, we chose the normalized transverse power of the field to be P = J \U\^d^ 
The numerical solutions are plotted in FIG. 2(a). 



cr = 2.5 
= 0.023. 



The power density of the soliton can be estimated by Ip — heoc 



h^ni 



Lp\p3l\ 



-P, where Lp is the transverse beam- width 

Moreover, the 



of the probe. Then the photon flux which passes through the soliton can be estimated by Ip/htop [1 

group velocity of the probe can also be estimated by [1^ : Vg « |^ , \ ^l^y 

we chose fin = 27, one can obtain that Ip = 1.81 x lO^"' mW/cm^, photon flux is 



They are all in proportion to \fl 



c\ 



If 



C = ^7, one can ootam tnat ip = i.oi x lu " mw/cm", pnoion nux is: 7 mm^ ns~ , and group velocity 
is: 3.15 mm/s, which have reach the nonclassical light region and the soliton can be termed as a NLS. 

If the detuning Ai > 0, Ep will act as self-defocusing in the medium. With the periodic modulation induced by 
fls, Eq. (10) has tt out-of-phase bright NLS (at the edge of the Brillouin zone) solutions. Such NLSs can be obtained 
by angle incident with the Bloch momentum lies in the vicinity of the edge of the flrst Brillouin zone. Diffraction in 
this condition is negative and self-defocusing nonlinearity is needed in trapping the flelds [y, [lOl- FIG. 2(b) plots a 
proflle of a TT out-of-phase DS with k = 2.5 (Ai = IOO7), Vb = 0.75 and P = 0.0106. 

In conclusion, we demonstrate that the optical induction technique can be applied to build an optical atomic 
discrete system in an atomic system through the technique of EIT. The probe field is experiencing the Kerr nonlinear 
self modulation and periodical index modulation of the lattice-forming wave, respectively. By changing the detuning 
Ai, both in-phase and tt out-of-phase NLSs can be produced in this system. The technique illustrated in this paper 



involves manipulating resonant optical properties of an atomic medium, and allow people to coherently control the 
propagation of a nonclassical light beam. By using this technique, many complex 2D LS including qusi-crystals LS 
[l7| . defect mediated LS [1^, honeycomb LS [l^, vortex i20!| or vortex-ring [2l| LS via optical induction can all be 
realized in such system. 

Yongyao Li thanks for Prof. Jianying Zhou for useful discussion. This work is supported by the NSFC (10934011, 
10774193, 10947126, ZA091201). 

Appendix A 

Here, we provide the detail deduction of the results in Eq. (4)- (7). According to the Eq. (3) and the initial 
condition: 



p(o) = 



/ 1 0\ 





\0 0/ 



(Al) 



the slowly varying amplitude evolution equations of the one-step density matrix elements are read as: 

pii^ = -^LD4lpil^ (A2) 

P42^ = -i^i2P^l2 - 2^cpi3\ (A3) 

pW = ~zCo,Jl^ 'ifhP^^l (A4) 

pi2 = -iJll (A5) 

4V = -»'^3iP^V + ^f^cP^i' + 'iflppf^ , (A6) 



P32 = -«'^32P32 " o^c(P33 " P22 )> (A7) 



P);^ = -73P^y + 2(^CP^3^ - ^hPW^ (A8) 

p'i^^-i^2j2i+\^*cPzi. (A9) 

p'i^ = -12P^^ + ^(^CP^S^ - ^CP^3^ (AlO) 

p'l^^-lJll (All) 

P^-pT^ (A12) 

where: 0)41 = 5 + A2 -4741, W31 = Ai -^731, W32 = A^- J732, ^^21 = 5 -1^21, and 7^^ = (7^+7^)72. Here, we assume 

71 w 72 ~ and 73 ~ 74 = 7. Because of C021 ~ 0, from Eq. (A6) and (A9), the steady state solution of p^i and P21 
are: 

(1) _ 2u}21^P _ *(1) 

^31 -4(:;3i^2i-|f^cP"^" ^ 

(1) _ »^»p _ ,(1) _ r^^»p 

^21 - 4^31^21 -If^cP^^^' ~ I^CP- ^ ^ 



And from this set of equations, one can easily obtain other steady state solution of one-step density matrix elements 
are all equals 0. Then the results in Eq. (4) can be obtained. 

The slowly varying amplitude evolution equations of the 2-step density matrix elements are as below: 

pfi' = -z(:;4ipii' + \^sP^2i. (A14) 

p§ = -iCoA2Pi2 - ipcP^iz, (A15) 

P43' = -i^A3P§ - if^cP^ih (A16) 

pfl = -lAP^il (A17) 

pf^ = -iCo^Ji^ + Ulcp^^l (A18) 

pf^ = -zCo.^p'i^ '-^c{pf^ pif) + '-p[\\ (A19) 

43 = -13P33 + \{^CP^2^ - ^*cpf^), (A20) 

p'i} = -iC02j2i + \^hp3i. (A21) 

P2^ = -12P^22 + \ i^*Cp32 - ^CP^23 ), ( A22) 

P[\^ = -12P[1\ (A23) 

f^-P^P- (A24) 
From Eq. (A20) and Eq. (A22), the level |3) and |2) with ilc can be treated as a two-level system. Then we have: 

M2 = P33' - P22 = -r32 {Wi^ - VKi°) ) + lincP^23 - ^*CP32 ) , (A25) 

where r32 « 73 is the transverse relaxation rate between |2) and |3) and W^2 ~ P33 ^ P22 = 0- Therefore, the steady 

(2) 
state solution of T432 can be obtained by: 

w(2) _ ^(2) (2) _ i .(-, (2) ^* (2). . ■ „„x 
i 32 

Substituted Eq. (A26) into Eq. (A19) and its conjugate, under the steady state condition, one can obtain an equations 

I W32P32 - a^^cH/ga - 2"PPl2 /Ar,7N 

t'^32P23 ~ 2^^C'^32 —2^^pP21 

(2) 

The steady state solution of pgj can be obtained from this equations set: 

(2) _ *(2) _ \np\^ncCo32 ,. . 



(2) (2) 

set about ^^32 and /?23 ^s below: 



From Eq. (A14) the steady state solution of p\^' can be obtained by: 

(2) _ *(2) _ il'^ilpfls 



Under the condition of Ai, A2 3> 732, |f^cp7 P32 ^^^ Pii '^^^ t)C simplified to: 

(2) _ \np\^nc (2) ^, n*(jnpns , , 

Other steady solutions of 2-step density matrix element are all equal to 0. Then the result of Eq. (5) can be obtained. 
Because we only needs the steady solution of pg-^ and p^2 ii^ the 3-step density matrix elements. Therefore in the 



3-step equations, we only need 2 equations sets: 

(3 
-IUJA2P42 

-iCj^Ji - i(f^&pi2^ + ^*ppf^ - ^sp'ih 



■~ „(3) jQ ^(3) 
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and 



(3) io (3) 
-«^3lP31 ~ 2"CP21 



>r/ = -^^^21^^!' + §(f^cP31^ - "PP23^ + f^sPfl^ 



(A32) 
By solving these two equations sets, one can obtain the results in Eq. (6) and Eq. (7). 
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